A phenomenology for multiple phases in the heavy fermion skutterudite 

superconductor PrOs4Sbi2 
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The two superconducting phases recently discovered in skutterudite PrOs4Sbi2 are discussed by 
using the phenomenological Ginzburg-Landau theory based on the the cubic (tetrahedral) crystal 
symmetry Th- Building on experimental input coming from the recent thermal conductivity mea- 
surements, the structure of the gap functions for these two phases are considered for the spin singlet 
pairing case. One of the phases, lying in the high temperature and high field region (A-phase) is 
a strongly "anisotropic s-wave" state, whereas the other phase (B- phase) is an "anisotropic s + id- 
wave" state. The B-phase breaks cubic crystal symmetry and time reversal symmetry spontaneously. 
The cubic symmetry breaking in B-phase causes anomalous crystal strains which could be seen in the 
anisotropy of thermal expansions and isothermal elastic constants. As a signal of the time reversal 
symmetry breaking, it is expected that an internal spontaneous magnetization should be observed 
by the /xSR measurement. 



I. INTRODUCTION 



Recently, heavy fermion superconductivity was discov- 
ered at 1.85 K in the skutterudite compound PrOs4Sbi2 
with a cubic (tetrahedral) crystal symmetry T/j [1]. In- 
triguingly, it was indicated that there are multiple super- 
conducting phases in this material, like in the Uranium 
compounds [2,3]. Two specific heat jumps were observed 
at Tci = 1.85i<r and Tc2 = 1.75i^ in the absence of the 
magnetic field [4]. Moreover, the thermal conductivity 
measurement in a magnetic field rotated relative to the 
crystal axes reveals a novel change in the symmetry of 
the gap function in the a6-plane [5]. The result was that 
the gap function has four point-node-like structure (four 
dips) in the [100] and [010] directions on the Fermi surface 
in the high field and high temperature region (A-phase), 
and two dips in the [010] direction in the low field and 
low temperature region (B-phase). The proposed phase 
diagram is given in Fig. I. Unfortunately, the gap struc- 
ture along the c-axis remains unresolved within in the 
presently available data [5], and one of the purpose of 
this paper is to consider this problem. The presence of 
the point node like structures appears to be consistent 
with the power law behavior of the specific heat [1,4] 
and the nuclear-spin-lattice relaxation [6] at low temper- 
atures (See also Ref. [7]). Another evidence for the dou- 
ble superconducting transition is the thermal expansion 
measurement along [100] direction [8]; two clear jumps 
are observed at around Td = l.^hK and Tc2 — 1.75K in 
the absence of the magnetic field (Fig. I). 

In this paper, the symmetry of the gap functions for 
A- and B-phase are discussed by use of a phenomeno- 
logical Ginzburg-Landau (GL) theory [2] based on the 
tetrahedral crystal symmetry T^. At present, there are 
no information as to whether the superconductivity is 
spin singlet or spin triplet pairing. Only the singlet pair- 



ing case will be discussed in this paper. The gap struc- 
ture along the c-axis is considered. We find that the gap 
function appears to have two additional dips along the 
c-axis in A-phase (in total there are six dips, along the 
[100], [010] and [001] directions) , while no additional 
dips are inferred for the B-phase resulting in a toal of 
two dips, along the [010] direction) . The A-phase gap is 
essentially an " anisotropic s-wave" state and the B-phase 
gap is an "anisotropic s + id- wave" state. The B-phase 
gap breaks the cubic crystal symmetry and time rever- 
sal symmetry spontaneously, from which one can expect 
a crystal deformation with cubic symmetry breaking to 
occur in this phase. To verify this, the phenomenological 
strain-order parameter coupling is taken into account and 
the thermal expansions and isothermal elastic constants 
are calculated [9]. One finds that these values become 
anisotropic in B-phase. 
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FIG. I. A rough sketch of the phase diagram proposed 
by Izawa et. al [5]. See, also Refs. [4] and [8]. 
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The pairing mechanism for this superconductivity 
should be also interesting. It is widely believed that the 
mechanism for the heavy fermion superconductivity in 
the Uranium compounds comes from magnetic interac- 
tions, but PrOs4Sbi2 is a non-magnetic compound. Sev- 
eral experiments indicate that the ground state of Pr 
4/-electrons in the crystalline electric field would most 
likely be the non-Kramers doublet T3 state [1,4,7]. The 
Ts ground state has quadrupolar fluctuations and its re- 
lations to the heavy fermion behavior [10] and also to a 
pairing mechanism seems to be possible to discuss. How- 
ever, its relation to the pairing symmetry, which will be 
focused in this paper, is unclear at present. Moreover, 
the singlet Fi state without quadrupolar fluctuations is 
also proposed as the ground state [11]. The details of the 
pairing mechanism will not be discussed here. 

The organization of this paper is as follows. In section 
II, the gap structures for the A- and B-phases are dis- 
cussed. In section III, jumps of the specific heat at two 
transition points (Td and Tc2) are calculated. To fit 
these calculations to the specific heat measurement [4], 
one obtains information for the coefficients of the fourth 
order terms in the GL potential. In section IV, argu- 
ments for a spontaneous cubic symmetry breaking in the 
B-phase are presented, and a summary is given in section 
V. 

II. THE GAP STRUCTURES 

In this section, the gap structures for the A- and 
B-phases in the spin singlet pairing case are discussed 
within the frameork of the phenomenological GL theory. 
Information on the gap structure along the c-axis, which 
is still unresolved by the thermal conductivity measure- 
ment [5], are obtained. 

The GL theory to be considered is based on the cubic 
crystal symmetry group Th- This group has essentially 
three irreducible representations Ai (Fi), E (F2,F3) and 
T (Ls). 

A. A-phase 

In this phase, four point-node-like (four-dip) structures 
of the gap function located within the a6-plane are ob- 
served along the [100] and [010] directions [5]. Then, the 
question is whether the gap function has two additional 
dips along c-axis (i.e. the gap has a total of six dips) or 
not (i.e. four dips in all). 

It should be pointed out that B.c2 is isotropic in this su- 
perconductor [5]. In general, the superconducting states 
represented by a multicomponent order parameter should 
have an anisotropic Hc2 [2,3]. Therefore, the gap func- 
tion in the A-phase appears to be represented by a single 
component order parameter, i.e., the Ai representation. 



The gap function with six dips along the [100], [010] 
and [001] directions belongs to the Ai representation be- 
cause this structure has the cubic symmetry T^- Actu- 
ally, one can see this as follows. Let rj and (^,,(k) stand 
for the order parameter and the basis function of the 
Ai representation, respectively. By the definition, (k) 
should be invariant under T^. The gap function in this 
representation is 

The six-dip state can be represented by using invariants 
[12]: 

M^) = s + kt + k^ + kt (s^-1), (1) 

and also 

M^) = s + klkl + k^kl + klkl (s«l), (2) 

• • • etc. are possible. When one takes s = — 1 in Eq. (1) 
(s = in Eq. (2)), the dips become point nodes. Gap 
functions such as in Eqs. (1) and (2) have been discussed 
in the context of the YNiB2C2 superconductor [13]. It 
should be noted that the function 0s (k) which can give 
rise to six dips is not limited to these two forms. There 
are many possible forms which have six dips. The point 
is that the six-dip state can be described in the Ai rep- 
resentation and essentially categorized into the so-called 
"anisotropic ,s-wavc" state. 

The four-dip state does not belong to the Ai represen- 
tation because this gap structure does not have the cubic 
symmetry. Actiially, one would see in the next subsection 
that this state belongs to the Ai©E combined represen- 
tation [12]. Then, a multicomponent orderparameter is 
needed to describe this gap fimction (See, Eq. (4)), and 
in general, an anisotropy in Hc2 is to be expected [2,3]. 
in contradiction to the experimental result [5] . Moreover, 
in the next subsection one would see that the four-dip 
structure arises in a somewhat accidental way (See, Eq. 
(5)). 

Then, in the A-phase, it seems to be natural to con- 
sider that the six-dip state (strongly anisotropic s-wave 
state) is realized. 

B. B-phase 

In this phase, two dips of the gap function within in 
afe-plane along [010] direction are observed [5]. Then, the 
question is whether the gap function has two additional 
dips along the c-axis (i.e. the gap has in total four dips) 
or not (i.e. a total of two dips). 

Phase transitions are observed by varying H and also 
varying T without a magnetic field turned on (Fig. I). 
To describe the phase transition without external fields, 
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one can consider the combined representations. Possi- 
ble combinations in the present case are Ai©E, Ai© T, 
EeT, and AiS E T. 

As we will see below, both the four-dip and two-dip 
states belong to the Ai® E representation. Let and 
5e stand for the effective coupling constant of the pairing 
interaction for the Ai channel and the E channel, respec- 
tively. In general, these values are split by the crystal 
field. The Ai® E combined state is realized when [2,14] 

9Ai ^ ffE, (3) 

is satisfied. This comes from some accidental degeneracy. 

The gap function in Ai® E model is expanded by the 
basis functions of Ai and E representations: 

A(k) = r](l)s{\s.) + m'Ply^-z-'-x^i^) + ^20z2_^2(k). (4) 

The function 4>2y'^-z'^-x'^(}^ is the first component of the 
2D basis of the E representation which transforms like 
the function 2ky — — k^. under the Th. group. It has 
two line nodes parallel to the ac-plane. The function 
(/)22_2,2(k) is the second component of the 2D basis of 
the E representation which transforms like the function 
fcf — fc^ under the T/j group. It has two line nodes crossing 
at two points along the [010] direction. Note that these 
functions may include terms beyond quadratic order in 
the momentum k. However the nodal structure of these 
basis functions remain unaltered, in general, even if the 
higher order terms are included. 

To obtain a four-dip state along the [010] and [001] 
directions, one must fine-tune the three basis functions 
into the particular forms 0s (k) = 1 — 2(A;^ + + 

(t>2y^-z^-x^{^) = 2fcy-fcf-fcf, <Az2_^2(k) = \/3(fc^-fc^), 

where the second order terms of k arc absent, and also 
the three order parameters are fixed to take the particu- 
lar values [12] 

(,y,r;r,r,2):.(l,-i ^). (5) 

This looks highly accidental. The gap function with four 
dips along [100] and [010] directions discussed in the pre- 
vious subsection can also be obtained in the same manner 
[12]. 

Two-dip state is realized when 

(?7,r?i,%) = (a,0,±i/3), a,/JeR, (6) 

because (t>^2_^2(}i) is zero and (/'s(k) becomes minimum 
along [010] direction. Note that the two-dip structure 
appears generically for all values of the two parameters 
a and (3. Moreover, it is not necessary to fix the basal 
functions 0s(k), (l)'2y^-z^-x'^{}^) and 0;j2_3.2(k) to partic- 
ular forms. The only assumption needed here is the six- 
dip structm-e of the function 0s (k) . It should be noted 
again that the function 0s (k) is not limited to a partic- 
ular form to obtain six dips. There are many possible 



forms to obtain six dips and described in the Ai rep- 
resentation (anisotropic s-wave state). Eq. (1) (or Eq. 
(2)) is one of the examples. Then, the state Eq. (6) is an 
anisotropic s + id^2_^2--weLve state and breaks the cubic 
crystal symmetry and also the time reversal symmetry 
spontaneously. 

Let us discuss the stability of these two states Eqs. (5) 
and (6). It will be shown that the state Eq. (6) is found 
as a stable state of the system, but the state Eq. (5) 
would be not. 

Consider the GL potential for Ai® E representation of 
T,. [14] 

J^op = Aa, {TM' + (3\vf + AE{T){\mf + 

+f3l{\Vl\' + \V2\'f + 02{VlV*2-V*lV2f 

+0i\v\'{\Vi\' + \V2\') 

+Y{^*\vi + v'2) + c.c}, (7) 

where 

^A,(T) = ^-1, Ae(T) = ^-1, 

and 

Ta, ~ exp[-l/5AiA^(0)], Te ~ exp[-l/5EAf(0)], 

in the weak coupling approximation (See. Eq. (3)). A^(0) 
is the density of states at the Fermi surface. We consider 
a temperature region T < mm[TAi,TE]. The terms pro- 
portional to 9i and 62 denote the coupling between the 
order paremeter in the Ai representation and the two 
order parameters in the E representation. In general, 
another coupling term 

could be included in the potential but this term makes 
the phase transition 1st order [14], which appears to be 

inconsistent with the specific heat experiments [1,4]. So, 
this term is excluded here. For simplicity, we take 

ei = 02. (8) 

By using the parameterization 

hi; Vi = \fi\ COS ipe^'^'- , 772 = I77I sinV'e"^^, 
Eq. (7) becomes 

Fop = AaAtm' + A^{TM' + PM* + 

-/?2|??|Sin2 2V'sin2(0i - ^2) + ^ihl'l??!' 
+^2|??n'?P [cos^ V cos 201 -I- sin^ V cos 2^2] • 

After a bit tedious calculation one obtains all of the sta- 
ble solutions [14]: 
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1. < /J2 



{v,VuV2) = ia,^e't±^^e^'^), (9) 



where, 



2(3 ' y 2(/3i- 1/321)' 
and is arbitral here. 
2. /32 < 

(c, ibirfcos^, ibirfsini/)) ,for < ^2 



where, 



{c,±d cos tp,±dsmip) ,for ^2 < 0, 

(10) 



c = 



Aa,(T)^^^ / ^e(T) 



2/3 



2/3i ' 



and ^ is arbitral. The degeneracy of tp is lifted by 
taking into account the 6-th order terms for (771,7/2) 
as a weak perturbation [2] 

7i(lml' + \v2\'f + 72{\vi\' + \m\'M + 4f 

+73|r/i|'|3r?2-»?il- 

The coupUng terms between 77 and (771,772) are ne- 
glected here. Then, 



0,7r (0 <73) 

±f (73 < 0). 



(11) 



From Eqs. (4), (6), (10) and (11), one can see that the 
two-dip state (anisotropic s + idz2_^2-'wave state) is sta- 
ble when /32,73 < < 62- On the other hand, it is hard 
to see how the four-dip state Eq. (5) can be stabilized in 
the present argument. 

The gap structures for the other stable solutions (Eq. 
(9) and Eq. (10) for 62 < 0) are inconsistent with 
the thermal conductivity measurement [5]. Furthermore, 
based on a consideration for the condensation energy 
within the weak coupling approach, it has been pointed 
out that the s + id-wave state is energetically favored in 
Ai (B E combined representation [15]. 

Fron these considerations, we conclude that the two- 
dip state appears to be realized in the B-phase. 



C. phase transition and symmetry breaking 

Let us summarize the results in this section and try to 
describe the phase transition in the absence of the mag- 
netic field. 



Put the parameters 

Ta, = Tel = 1.85K, Te = Te2 = 1.75K, (12) 

and /32,73 < < ^2 in the GL free energy Eq. (7). For 

simplicity, wc take Eq. (8). Then, the stationary values 
of the order parameters in the A-phase {Tc2 <T < T^) 
are 



^A,(r) „ 



(13) 



and in the B phase (T < Tc2) 



0,7/2 = ±M 



^e(T) 



2/3i 



(14) 



From Eqs. (4), (13) and (14), one can see that the gap 
function in the A-phase has six dips along the [100], [010] 
and [001] directions (anisotropic s-wave state), and that 
in the B-phase there are two dips along the [010] direc- 
tion (anisotropic s + idz2_x2-wa,ye state). The result is 
consistent with the thermal conductivity measurement 
[5]. 

The gap function in the B-phase breaks the cubic crys- 
tal symmetry and also the time reversal symmetry. The 
consequence of this cubic symmetry breaking will be dis- 
cussed in section IV. As a signal of the time reversal 
symmetry breaking, it is expected that a spontaneous 
magnetization be observed by the fiSR measurement [2]. 



III. SPECIFIC HEAT JUMP 

In this section, we calculate the jumps of the spe- 
cific heat at the two transition points Td = 1.85K and 
Tc2 = 1.75K in Fig. I indicated by the specific heat mea- 
surement [4]. To fit the calculation to the experimental 
result the ratio of the coefficients of 4-th order terms in 
GL potential /3//3i is obtained. The simplification Eq. 
(8) is used. 

Let us introduce 

V= iv,Vi,V2)- 

By using the stationary condition for the free energy Eq. 
(7) 



df] ' 



(15) 



the specific heat is 



C(T) = Co{T) - T 



dT2 



„ ^ r Or] d^F 

^?o(T)-T|— 



(16) 
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where Co{T) is the background specific heat. Let us de- 
fine AC(T) = C(r) - Co(T). Then, in the A-phase, 



AaC(T) 



while in the B-phase, 



AbC{T) = T{,^+ ^ 



2/?r2 2/31T3V ■ 



(17) 



(18) 



Therefore, the ratio o£ the jumps at the two transition 
points is 



ABC(Te2) - AAC(Te2) _ /JT^i 



(19) 



The experimental result for this quantity is around 1 [4]. 
Then, from Eq. (12), one can obtain the ratio 

13/ Pi ~ 0.9. (20) 



IV. CUBIC SYMMETRY BREAKING IN 
B-PHASE 

the B-phase gap breaks the cubic crystal symmetry Th 
spontaneously. So, one can expect a crystal deformation 
with the cubic symmetry breaking. In this section, the 
strain-order parameter coupling is taken into account and 
thermal expansions and isothermal clastic constants are 
calculated [9] . In the B-phase, one obtains anisotropic re- 
sults for these values. This comes from the spontaneous 
crystal symmetry breaking. 



A. strain-orderparameter coupling 

The GL potential which include the strain-order pa- 
rameter coupling is [2], 



Fql = ^op + -Fst-op + -Fell 

-Fst-Op = -\'nfCAi{<^XX + eyy + Czz) 

-|?7l|^ [CA^{exx + eyy + izz) 

+ CE(2eyy — Czz — ixx)] 

-M'^ [CAi (Exx + £3/3/ + ezz) 



CE(2e, 



)] 



-•s/3Ce(?7i?72 -I- c.c.){ezz - e^x) 
-C'aeIwi + c.c.){2eyy - Czz - e^x) 

-V3Cae(W2 + C.C.){ezz - Exx), 

^ei = ^Ei°i/^'i^«' (21) 

ijkl 

where c^°],; is the background elastic constants. Obvi- 
ously, it is symmetric under the Th symmetry. 



Let us consider the amplitude of the strain. It is ob- 
tained by solving the equation 







dF, 



GL 



We introduce the notation e^x = ^i^^yy = £2,£zz = 
esjExj/ = e4,€yz = €5,ezx = £6- In the A-phase (see, Eq. 
(13)), one finds the isotropic result 



ei = £2 = £3 = - 



64 = £5 = £6 = 0, 



A 2/3 '■ 



(22) 



where A = c^^ = C22 = 0^33 ■ In B-phase (See, Eq. (14)), 
an anisotropy can be seen in the y direction: 



£3 = £1 



Ca.AaAT) Ca,+CeAe{T) 
' X 2(3 A 2/3i ' 

Ca.AaAT) Ca,-2CeAe{T) 



'' = -X 2(3 

£4 = £5 = £6 = 0. 



2/3i 



(23) 



It should be noted that the terms proportional to (7ae 
in Eq. (21) play no roles in the state Eq. (14). 



B. thermal expansion 

The presence of the strain can be measured by the 
thermal expansion which is definition by 



where ji = 1,- • -,6, and aj}''' is the background value. 
Define Aa^ = a^ — a|f ' . Then, in the A-phase (See, Eq. 
(22)), 

AaOi = AAa2 = AaQ;3 = A, 

A /pi cl 

AAa4 = AaQs = AA^e = 0. (25) 

In B-phasc (Sec, Eq. (23)), the thermal expansion be- 
comes anisotropic: 



(24) 



Aeas = Aeai 

ABa2 : 



Cai 1 Cai -I- Ce 1 



A 2(3Tei A 2(3iTc2 

Cai 1 Cai — 2Ce 1 



A 2(3Tci A 2(3xTo2 

ABa4 = Aeas = Aeae = 0. (26) 

The thermal expansion along [100] direction has been 
measured [8], and two clear jumps are observed at Td 
and Tc2. The ratio of these two jumps is 
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From Eqs. (12) and (20), one can obtain a relation 

Ce/Ca, ^ 0.9, (28) 

and estimate the ratio of jumps at the two transition 
points for the thermal expansion in the [010] direction 

^H^P^ = fl-2-^U^.-0.8. (29) 

This quantity has not been measured yet. 



C. isothermal elastic constants 

The isothermal elastic constant is 
'd-'FGL\ 



(30) 



The small variation of the strain changes the amphtude 
of the order parameters as 



6\v\^ = ^{Se,+6e2 + 6es), 



-^(2<5e2-<5e3-<5ei). 
By using the condition Eq. (15), 



c = c(°) 



define Ac^i, = c^i, — c^^},. In A-phase, 



c.c. 



(31) 



(32) 



where c|S is the background clastic constant. Let us 



(J2 

AACn=AAC22 = --- = AAC31 = -^, (33) 

and the others are zero. In B-phase, 

Cl^ (Ca,+Ce)^ 



AbC33 = AbCh = AbCsi 



2/3 



2/3i 



A ^ ^Ai (Cai - 2Ce)^ 

AbC22 = -^ ' 



(34) 



AbC23 = AbCi2 = - 



Cl^ (Cai+Ce)(Cai-2Ce) 



2/3 



2/3i 



and the others are zero. From Eqs. (20) and (28), one 
can estimate the ratios of jumps at two transition points 
of the elastic constants, 



AbCh - AaCh 




/ 


Ce a 


Aach 


Pi 




AbC22 - AaC22 


/3 




2^ 


AaC22 


/3i 




Cai , 


AbCi2 - AaCi2 


/3 




Ce\ 


AaC12 


/3i 


CaJ 



77^ ^3.6, 



~ 0.6, 



,Ce_ 

CAi 



-1.4. 



(35) 



In conclusion for this section, the thermal expansion 
and the elastic constant are isotropic in the A-phase, but 
anisotropic in the B-phase. To observe the anisotropy 
would serve as additional evidences for the phase tran- 
sition and the spontaneous cubic symmetry breaking in 
the B-phase. 



V. SUMMARY AND DISCUSSIONS 

In summary, the multiple superconducting phases (A- 
phase and B-phase) in the skutterudite PrOs4Sbi2 with 
cubic crystal symmetry T;i is discussed by using the phe- 
nomenological Ginzburg-Landau approach. The spin sin- 
glet pairing case is considered here. According to the 
thermal conductivity measurement [5] , the gap function 
has within in the a6-planc four point-node-like structure 
(four dips) along the [100] and [010] directions on the 
Fermi surface in the A-phase, and two dips along the 
[010] direction in the B-phase. But the gap structure 
along the c-axis was unresolved. This problem is consid- 
ered here. 

Because of the absence of an anisotropy in IIc2 [5], 
the superconductivity in A-phase appears to be the con- 
ventional state (anisotropic s-wave state). Then, the gap 
structure has the cubic symmetry (the Ai representation) 
and leads to two additional dips along the c-axis. In the 
B-phase, a state without additional dips along the c-axis 
(anisotropic s + «rf^2_2,2-wave state) is stabilized by the 
Ginzburg-Landau potential of the Ai®E combined rep- 
resentation in a wide parameter region. On the other 
hand, it is hard to see that a state with two dips along 
c-axis is stabilized, since to obtain such a state requires 
a highly accidental mixing between the basis functions of 
the crystal symmetry T/j. Moreover, it has been shown 
that the s + id-wave state is energetically favored in the 
Ai©E combined representation based on the considera- 
tion for the condensation energy of each state within the 
weak coupling approach [15]. Hence it is natural to ex- 
pect that the gap function in A-phase has a total of six 
dips along [100], [010] and [001] directions and while in 
the B-phase there are in total two dips along the [010] 
direction. The phase transition can be described by the 
Ginzburg-Landau potential. 

The gap function inthe B-phase breaks the cubic crys- 
tal symmetry T/j spontaneously. This symmetry break- 
ing would cause anisotropic crystal strain and anisotropic 
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anomaly of the thermal expansion and the isothermal 
elastic constant could be observed . 

The gap function in B-phase breaks the time reversal 
symmetry (See, Eqs. (4) and (6)). Then, the sponta- 
neous magnetization around impurities is expected to be 
observed by the /xSR measurement [2]. The magnetiza- 
tion have not been detected yet. More accurate measure- 
ments are highly desired [7]. 

Another expected phenomenon in the B-phase is 
anomalous flux flow which comes from the domain for- 
mation [16], because the domain structure would exist 
in the B-phase since the stable state Eq. (6) is two-fold 
degenerate. 

Additional phase transitions induced by uniaxial pres- 
sure may be possible in this superconductor. Such a tran- 
sition could occur in a superconductor described by mul- 
ticomponent order parameter [2]. 
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